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(T) Progress in computing the hadron spectrum in lattice QCD using stochastic LapH quark propaga- 

tors is described. The stochastic LapH algorithm is a particular quark smearing algorithm that also 
£> allows the computation of all-to-all quark propagators. All-to-all quark propagators are required 

in our approach of using a large set of spatially extended hadron operators and explicit multi- 



5_i particle operators to access excited states. We report on the progress made in the various isospin 



channels on 2+1 dynamical, anisotropic lattices generated by the Hadron Spectrum Collaboration. 
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1. Introduction 

One of the goals of our collaboration has been to compute, from first principles lattice QCD 
simulations, the low-lying hadron spectrum in various symmetry channels in hopes to aid our un- 
derstanding of QCD through comparisons with low energy effective theories. The reliable extrac- 
tion of the energies of the excited states in each symmetry channel requires that we use interpolating 
operators which couple well with the ground states in the channel and also reasonably well with the 
excited states up to the energy scale that one is interested in. The progress made in the last couple 
of years for single-particle and some two-particle operators/correlation functions have appeared in 
various papers and proceedings ([l]-[7]). 

The extraction of excited energy levels on dynamical configurations must confront the problem 
of mixing of single-particle states with multi-particle states and requires new techniques to be 
developed given the finite resources available. The construction of a set of good interpolating 
operators which have good signal-to-noise ratio is crucial in computing the excited states through 
the variational method as the signal for the heavier excited states decay rapidly into the noise, 
destabilizing the diagonalization procedure. In this contribution, we report on our progress towards 
solving these issues. 

2. Lattices 

The dynamical lattices that we use are the clover-improved, anisotropic lattices ([8]) that have 
been generated by the HadSpec collaboration [9]. The anisotropy not only allows one to extract 
energies at early timeslices (as long as the interpolating operators are effective) before the signal- 
to-noise ratio deteriorates but also generally has smaller fluctuations. There has only been a single 
lattice spacing generated so far corresponding to ro/a s = 3.221(25). The anisotropy was tuned to 
3.5 with a measured valued of 3.3 using the pion dispersion relation. The quark mass used in this 
study corresponds to pion masses of 240 MeV (250 configurations). 

3. Interpolating Operators 

The interpolating operators that are used to form the correlation matrix for digaonalization is 
constructed so that they transform irreducibly under the lattice symmetry group of rotations and 
translations and have the appropriate flavor structure and G-parity. Operators for hadrons at rest 
can be constructed from a handful of conventional quark propagators as long as the displacement 
types are restricted to the simple displacements. The interpolating operators that we will be using 
are based on the irreducible representations of the appropriate cubic point groups and the pruning 
of these types of operators have been studied in detail in previous studies. The choice of the 
elemental operators to use were based on its simplicity and diversity while maintaining the best 
overlap and low noise-levels in the single particle correlation functions. A preliminary study of the 
isovector meson spectrum was presented last year at the lattice conference for the heavier quark 
mass corresponding to m % = 390 MeV. 

The same approach is taken in designing/pruning the interpolating operators for hadrons with 
non-zero momenta, k (in units of 2n/L). In this study, we consider those states which have mo- 
menta of the form, k = (0,0, n), k = (0,n,n) and k = (n,n,n). The group theory for these choices 
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of momenta was discussed in the conference report from last year [1] and the general discussion of 
the group theory is given in [10]. 

3.1 Quark Propagators 

In order to incorporate hadron operators with non-zero momentum and explicit two-particle 
operators in our operator set, we construct the hadron interpolating operators from stochastic LapH 
quark propagators ([1 1]). The stochastic LapH algorithm is an efficient all-to-all extension of the 
distillation algorithm ([12]) which does not suffer from poor volume scaling behaviour. The method 
allows one to inject any momentum into the interpolating operators after the propagators have been 
generated, the construction of extended operators and it also separates the source operator from 
the sink operator which greatly simplifies the construction of large conelation matrices needed to 
extract the excited state spectrum. 



3.1.1 Stochastic LapH Quark Propagators 

In this section, we briefly review the stochastic LapH algorithm. The source vector (for quark 
label A) is given by 

pgg'(x,0=(A? ) V,pMpW(«,0) v 

where b is the color index, j8 is the spin index, Dr*' is the covariant derivative in the k direction with 
separation p, V5 is the matrix of eigenvectors of the Laplacian and is the projection operator 
for the dilution scheme being applied. The dilution schemes can be categorized into four different 
types: 

V [ fJ = 8ij, d = 0, (no dilution) 

Pp = 8ij 8 di , d = 0,. . . ,N — 1 (full dilution) 
P,7 = Stj S dt [Ki/N\ d = 0,...,K-l, (block-*) 
pH ] = 5ij 8 di i mod K d = 0, . . . , K - 1 , (interlace-*) 

where N is the dimension of the space of dilution type and N/K is an integer. 

The solution vector, ^> a ajf, is then computed by inverting CI = 74M with this source. The 
solution vector is then smeared with the smearing operator and displaced appropriately, 

Vaajp ~ & V ^ ) aaM PWP 
where Q. = 74M and M is the Dirac operator. The quark propagator (for line A) is then given by 

d 

3.2 Meson Correlation Functions 

Meson correlation functions are formed by correlating meson operators given by 



(p 1 ,<p 2 ;p,0 = c^£^'>( x+ ( d ^))p^ 1 (x,?;p0*<p5 ] (x, ? ;p 2 ) 
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Figure 1: Preliminary results for the glueball corre- Figure 2: Preliminary results for the correlation 
lation function on the 24 3 lattice at 240 MeV pion function of the glueball and %-% mixing term for the 



mass. 



240 MeV configurations. 



where d a are introduced to maintain G-parity when needed. The general form of the meson corre- 
lator takes the form (many of the indices have been suppressed in this expression for brevity), 



Cu(p,t f -t ) 



SAA 8 BB^i dld2] {$h9i\V,t. 



+8 AB 8 As ^ d,d,] (pi , <pi ;p, 



(3-1) 



[d 2 d 2 ] 



((p2,p2-,V,t )* 



where the second term appears for isoscalar mesons. 

4. Preliminary Results 

In this section, we present a few selected preliminary results for particles in different isospin 
sectors. 



4.1 7 = 

The isoscalar channel is one of the most challenging channels as there are disconnected 
("same-time") diagrams and "box" diagrams appearing in the correlation function. The com- 
pletely disconnected diagrams generally suffer from large gauge noise and hence is an ideal case 
for stochastic estimation over 'non-stochastic' evaluations such as the distillation algorithm. The 
"box diagrams" are just as difficult to evaluate but are perhaps not as noisy. The stochastic LapH 
algorithm allows the computation of all of these diagrams which are involved in this channel from 
four independent noise sources. Two noise sources are needed just to construct a single pion opera- 
tor and so at the expense of computing two more noise vectors, we can evaluate all of the diagrams 
needed for the isoscalar channel. The scalar glueball operator comes for free as well because one 
can construct it from the eigenvalues of the LapH operator. A glueball correlation function using 
these LapH operators is shown in Fig. 2 on the 24 3 x 256 lattice with the 240 MeV pion mass. 
We have studied the mixings between the scalar meson and the scalar glueball as well as scalar 
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Figure 3: Preliminary results for the various contributions to the nn — fo mixing term. 



mesons with the 1 = kk scattering state. The off-diagonal correlation functions for these differ- 
ent interpolating operators for the 1 = channel on the 24 3 lattices are shown in Fig. 2 and Fig. 3. 



There are no disconnected diagrams nor any box diagrams that needs to be evaluated for kaons. 
The gain from the stochastic LapH algorithm is that different momenta values for the kaons can be 
simulated by introducing simple exponential factors into the zero momentum operators. However, 
there will be an extra noise source/solution pair needed at the strange quark mass. Our preliminary 
results suggest that the introduction of the stochastic noise solely in the LapH subspace does not 
make the determination of energies with finite momentum any harder than zero momentum states 
with similar energies. We show an example for the kaon with momentum (1, 1, 1) in lattice units in 



The rho decay width has attracted attention in recent years as several groups have performed 
large scale simulations at light quark masses ([13]-[15]). The lattice sizes and the quark masses 
we have are limited as larger lattices with lighter masses are still being generated now. We have 
generated many of the operators with finite momenta that will be needed for the analysis of this 
channel. We show one of the representative effective masses in this channel in Fig. 5. 



The isospin-2 channel is an exotic channel and does not require any disconnected diagrams 
to be evaluated. The correlation function can be evaluated by computing the product of two pion 
propagators and the corresponding quark-exchanged diagrams. Since the finite momentum single 
pion operators have already been computed, we can use many different Tin operators to reduce 
the contamination from excited states. The use of many operators for this state is important in 
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Figure 4: The effective mass of the kaon with momentum (1,1,1) in lattice units. 



a large volume because the level spacings decrease with increasing volume. We construct a cor- 
relation matrix which is then diagonalized to find the optimized correlation function. We fit the 
optimal correlator for various t m i n values to check for stability in the fits. A "?„„■„ plot" is shown in 
Fig. 6 after diagonalization. The scattering length that is extracted from the fits are consistent with 
continuum two-loop chiral perturbation theory predictions [16], but the value is sensitive to the 
choice of the optimization timeslice. We are examining the systematic errors of choosing different 
diagonalization times and also the sensitivity to t max included in the fits. 

5. Summary 

The stochastic LapH quark smearing algorithm has been tested for single particle meson states 
with various momenta and isoscalar channels including the mixing of the %-% state with the scalar 
glueball and the scalar meson state. The elastic scattering length was estimated using Liischer's 
formula and found to be consistent with large errors with continuum chiral perturbation theory. 
The 1 = 2 exotic state was computed using various finite momenta pion operators and a result 
consistent with chiral perturbation theory was obtained. The systematic errors due to the different 
diagonalization procedure has not been fully understood yet and is a subject of current investigation. 
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